Interaction effects on thermal transport in quantum wires 



O 

Or 



G 

O 

o 



> 
in 

cn 

On 
O 



X 



Alex Levchenko, 1,2 Tobias Micklitz, 3 Zoran Ristivojevic, 4 and K. A. Matveev 1 

Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA 
2 Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA 
3 Dahlem Center for Complex Quantum Systems and Institut fur Theoretische Physik, 
Freie Universitat Berlin, 14195 Berlin, Germany 
* Laboratoire de Physique Theorique-CNRS, Ecole Normale Superieure, 24 rue Lhomond, 75005 Paris, France 

(Dated: September 16, 2011) 

We develop a theory of thermal transport of weakly interacting electrons in quantum wires. 
Unlike higher-dimensional systems, a one-dimensional electron gas requires three-particle collisions 
for energy relaxation. The fastest relaxation is provided by the intrabranch scattering of comoving 
electrons which establishes a partially equilibrated form of the distribution function. The thermal 
conductance is governed by the slower interbranch processes which enable energy exchange between 
counterpropagating particles. We derive an analytic expression for the thermal conductance of 
interacting electrons valid for arbitrary relation between the wire length and electron thermalization 
length. We find that in sufficiently long wires the interaction-induced correction to the thermal 
conductance saturates to an interaction-independent value. 

PACS numbers: 72.10.-d, 71.10.Pm, 72.15.Lh 



I. INTRODUCTION 

The classical Drude theory of electronic transport pro- 
vides universal relation between electric and thermal 
transport coefficients known as the Wiedemann-Franz 
lawi 
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Here K and G are, respectively, the thermal and elec- 
tric conductances, T is the temperature in energy units 
(fee = 1), and e is the electron charge. The relation 
represented by Eq. QJ is a natural one since for non- 
interacting particles both charge and energy are carried 
by the electronic excitations. Furthermore, as long as 
elastic collisions govern the transport, the validity of 
the Wiedemann-Franz law has been confirmed in the 
case of arbitrary impurity scattering? 2 - However, an ac- 
count of the electron-electron interaction effects within 
the Fermi-liquid theory gives corrections to both G and 
K££ These lead to a deviation from the Wiedemann- 
Franz law [Eq. (JTJl], which is associated with the inelastic 
forward scattering of electrons. In general, violation of 
the Wiedemann-Franz law is a hallmark of electron in- 
teraction effects and thus is of conceptual interest. 

In one-dimensional conductors, such as quantum wires 
or quantum Hall edge states, the electron system can 
no longer be described as a Fermi liquid but instead is 
expected to form a Luttinger liquid. 5 It has been shown 
that for the perfect Luttinger liquid conductor, such as an 
impurity-free single-channel quantum wire, interactions 
inside the wire neither affect conductance quantization 6 



nor change the thermal conductance of the system^ 
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Since both G and K remain the same as in the case of 
noninteracting electrons, the Wiedemann-Franz law ([1]) 
holds for an ideal Luttinger liquid conductor. 

There are two important exceptions known in the lit- 
erature. The first one is a Luttinger liquid with an impu- 
rity studied by Kane and Fisher £ In that case electron 
backscattering takes place which strongly renormalizes 
both G and K, such that the Wiedemann-Franz law is 
violated. The second case is the Luttinger liquid with 
long-range inhomogeneities studied by Fazio et alZ If the 
spatial variations related to these inhomogeneities occur 
on a length scale much larger than the Fermi wavelength, 
electrons will not suffer any backscattering. The electric 
conductance will, therefore, be given by its noninteract- 
ing value [Eq. (J5J)]. At the same time the thermal con- 
ductance K will be altered by interactions. The reason 
for this is as follows. For the system with broken trans- 
lational invariance, momentum is not conserved. As a 
result, there are allowed certain pair collisions which con- 
serve the number of right and left movers independently 
but provide energy exchange between them. These are 
precisely the scattering processes that thermalize elec- 
trons and thus lead to the violation of the Wiedemann- 
Franz law. 

Recent advances in the fabrication of tunable con- 
strictions in high-mobility two-dimensional electron gases 
have allowed precise and sensitive thermal measurements 
in clean one-dimensional systems. These include ex- 
periments on the thermal transport of single-channel 
quantum wires, where a lower value of the thermal 
conductance than that predicted by the Wiedemann- 
Franz law was observed at the plateau of the elec- 
trical conductance^! 1 ^ Another set of experiments re- 



2 



ported enhanced thermopower in low-density quantum 
wire o 11 ! 12 and quantum Hall edges Remarkable ex- 
periments based on momentum-resolved tunneling spec- 
troscopy provided direct evidence for the electronic ther- 
malization in one-dimensional systems*""— Clearly in- 
teraction effects are responsible for the observed features; 
however, Luttinger liquid theory does not provide an ad- 
equate description for these observations. 

Ongoing theoretical efforts in the study of one- 
dimensional electron systems focus on nonequilibrium 
dynamic o 17 ' 18 and consequences of the nonlinear disper- 
sion in transport^— that lie beyond the scope of con- 
ventional Luttinger liquid paradigm. The kinetics of one- 
dimensional electrons with nonlinear dispersion is pecu- 
liar. Indeed, constraints imposed by the momentum and 
energy conservation allow either zero-momentum trans- 
fer or exchange of momenta for pair collisions. Neither 
process changes the electronic distribution function, and 
as a result they have no effect on transport coefficients 
and relaxation. Therefore, three-particle collisions play 
a central role^ 

In this paper we study the fate of the Wiedemann- 
Franz law and the origin of the electron energy relaxation 
in clean single-channel quantum wires, accounting for the 
scattering processes that involve three-particle collisions. 
The paper is organized as follows. In the next section, 
Secim we place our work into the context of recent stud- 
ies on equilibration in quantum wires and explain the 
concept of partially equilibrated electron liquids, which 
is central for our study. We then develop in Sec. IIIII 
the theory of the thermal transport in one-dimensional 
electron liquids based on the Boltzmann equation with 
three-particle collisions included via the corresponding 
collision integral and scattering rate. We elucidate the 
scattering processes involved, discuss the role of the spin 
and limitations of our theory, and summarize our findings 
in Sec. lIVI Additional comments and directions for future 
work are presented in Sec. |Vj Various technical aspects 
of our calculations are relegated to Appendixes (TA"l)- ([Cl) . 

II. PARTIALLY EQUILIBRATED 
ONE-DIMENSIONAL ELECTRONS 

A. Noninteracting electrons 





FIG. 1: (Color online) Schematic picture of a quantum wire 
of length L. Electrons in the left and right leads are described 
by Fermi distribution functions characterized by temperatures 
Tu T \ and chemical potentials iiu r y Due to three-particle col- 
lisions electrons may backscatter and also exchange energy 
between the subsystems of warmer right movers and colder 
left movers. 



T + AT/2, T r — T — AT/2, and m = fi + eV/2, 
fi r = ji — eV/2 are the different temperatures and chem- 
ical potentials of left and right leads, respectively (see 
Fig. [T]) . Employing the distribution function from Eq. (U) 
to linear order in V and AT, one readily finds the electric 
and heat currents / = GqV\at=o and Iq = KqAT\j = q, 
with conductances Go and Ko, which coincide with the 
earlier stated noninteracting values, Eqs. © and 27 

In the presence of weak interactions the distribution 
((3]) describes an out-of-equilibrium situation. Collisions 
lead to electrons exchanging energy and momentum, with 
some particles experiencing backscattering. As a result, 
net particle (N R ) and heat (Q R ) currents flow between 
the subsystems of right- and left-moving electrons, relax- 
ing V and AT (see Fig. [1] for a schematic illustration) . 
The effect of electron-electron collisions on the distribu- 
tion function depends strongly on the length of the wire. 
Short wires are traversed by the electrons relatively fast, 
leaving interactions only little time to change the distri- 
bution in Eq. Q considerably. In the limit of a very long 
wire, on the other hand, one should expect full equilibra- 
tion of left- and right-moving electrons into a single dis- 
tribution, even in the case of weak interactions. It turns 
out that there exists a hierarchy of three-particle scatter- 
ing processes, classified by the corresponding relaxation 
rates or, equivalently, inelastic scattering lengths, which 
have different effects on the electron distribution function 
Eq. Q. 



Noninteracting electrons propagate ballistically 
through the wire. They keep a memory of the lead 
they originated from and remain in equilibrium with 
the corresponding lead electrons. If voltage (V) and/or 
temperature differences (AT) are applied across the 
wire, right- and left-moving particles are at different 
equilibria, and the distribution function takes the form 

6(jp) 8(-p) 



f P = 



(4) 



e T i +1 e ?v + 1 

where e p = p 2 /2m is the energy of an electron with 
momentum p and 9(p) is the unit step function. T; = 



B. Partially equilibrated electrons 

We start our discussion of the different scattering pro- 
cesses involved in the electronic relaxation with the pro- 
cess shown in Fig. HJa). This three-particle collision 
provides z??ira- branch relaxation within the subsystem 
of right-moving electrons. A similar relaxation process 
for the left movers is not shown in the figure but is im- 
plicit. This process can be described by the correspond- 
ing inelastic scattering length, which we denote in the 
following as t a . The precise form of £ a and its temper- 
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ature dependence are model specific. They are deter- 
mined by the scattering amplitude for the given inter- 
action potential and the phase space available for this 
scattering [Fig. HJa)] to occur. Quite generally one may 
argue that l a scales as a power of T. Indeed, at low tem- 
peratures, T <C fJ>, all participating scattering states are 
located within the energy strip ~ T near the Fermi level. 
What is important for the present discussion is that for 
wires with length L ^> l a intrabranch electron collisions 
[Fig. [H[a)] become so efficient that the initial distribu- 
tion function Eq. (Q| will be modified by interactions. 
One can find the resulting distribution by employing the 
following observation. Intrabranch collisions conserve 
independently six quantities. These are the number of 
right and left movers, N R t L — E P >o fpi their momenta, 
P R/L = Ep ?0 P/p» and energies, V/ L = £^ e p / p . 



The form of the resulting electron distribution function 
f p can be obtained from a general statistical mechan- 
ics argument by maximizing the entropy of electrons, 
S = -E p [/p ln /p + (!- /p) m (! ~ f P )], under the con- 
straint of the conserved quantities^ 



fp = 



9(p) 



+ 



0(-p) 



(5) 
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This distribution is characterized by six unknown param- 
eters (Lagrange multipliers) which have a transparent 
physical interpretation. Indeed, in Eq. ([5]) T R / L are ef- 
fective electron temperatures for right and left movers 
different from those in the leads. The parameters u R / L 
have dimension of velocity and account for the conserva- 
tion of momentum in electron collisions. Finally, \i L l R 
are uncquilibratcd chemical potentials of left- and right- 
moving particles. In principle, all these parameters may 
depend on the position along the wire. 

For longer wires interbranch three-particle collisions 
[see Fig. EJb)] become progressively more important. 
Unlike intrabranch relaxation these processes allow en- 
ergy and momentum exchange between the subsystems 
of right and left movers; thus P R / L and E R I L are no 
longer independently conserved. However, the full mo- 
mentum P = P R + P L and energy E — E R + E L are 
obviously conserved. Corresponding to Fig. [2jb) , the 
scattering length ij, is model specific and calculated in 
Appendix [C] We note here that for all interaction po- 
tentials we studied lb/l a ~ fJ>/T >• I. In view of this 
distinct length scale separation, lb ^> l a , the electron 
distribution in Eq. <j5j> is established at the first stage of 
the thermalization process. However, for longer wires, 
L 3> lb, relaxation of counterpropagating electrons be- 
comes so efficient that the temperatures T R ^ L and boost 
velocities u R l L of right and left movers become equal, 
T R = T L = T and u R = u L — u, due to energy and 
momentum exchange. At the same time, the chemical 
potentials [i R l L are still unequilibrated in this regime, 
Afi = n R — [i L 7^ 0, since the numbers of right- and left- 
moving electrons are still independently conserved. As 
a result, for wires with length L 3> lb the distribution 





FIG. 2: (Color online) (a) Intra-branch relaxation of co- 
moving electrons that establishes the partially equilibrated 
form of the distribution function in Eq. <(Sj . (b) Dominant in- 
terbranch three-particle process for the energy exchange Q R 
between counterpropagating electrons that contributes to the 
thermal conductance correction, (c) Leading three-particle 
collision which results in a finite rate N R oc e _M//T and thus 
a temperature-dependent correction to the conductance of a 
short wire (Ref. [l9h . (d) Equilibration mechanism: multistep 
diffusion through the bottom of the band of an electron from 
the right to the left Fermi point accompanied by the excita- 
tion of many electron-hole pairs (Refs. l2CH22h . 



function ([5]) transforms into 
6{ P ) 



f P 
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Because A/i / Owe refer to the states of electron system 
described by the distributions in Eqs. (J5J) and © as the 
states of partial equilibration. 



C. Fully equilibrated electrons 

Energy and momentum conservation allow for the scat- 
tering process in which an electron at the bottom of the 
band is backscattered by two other particles near the 
Fermi level [see Fig.^c)]. This is the basic three-particle 
process that changes the numbers of right and left movers 
before and after collision. In particular, the exponentially 
small discontinuity of the distributions Eqs. ([5]) and © 
at p = will be smeared by collisions of this type. 

Complete equilibration of electrons, namely, relaxation 
of ApL, relies on the electron backscattering from the right 
to the left Fermi point. One should notice here that it 
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is impossible to realize such scattering directly since it 
requires momentum transfer of 2pp while Fermi blocking 
restricts typical momentum exchange in the collision to 
bp ~ T/vp <C pf- As a consequence, complete electron 
backscattering, and thus relaxation of the chemical po- 
tential difference A /it, occurs via a large number of small 
steps Sp in momentum space such that Sp <C Pf^ In its 
passage between the subsystems of right and left movers 
the backscattered electron has to pass the bottleneck of 
occupied states at the bottom of the band [see Fig.^d)]. 
As a result, backscattering of electrons is exponentially 
suppressed by the probability ~ e~ M / T to find an un- 
occupied state at the bottom of the band, and thus the 
equilibration length £ eq for the relaxation of the difference 
in chemical potentials of left and right movers is expo- 
nentially large, £ eq oc e^/ T . For sufficiently long wires, 
L l e q, the state of full equilibration is achieved and 
described by the distributio n 2 ^ 21 

fp = e p -v d p-»e q ' CO 

e 1 + 1 

where the chemical potential /j, eq inside the equilibrated 
wire is, in general, different from /i/( r ) in the leads. In 
Eq. Vd = I/ne is the electron drift velocity, where 
I is the electric current and n is the electron density. 
The partially equilibrated distribution function given by 
Eq. ((6l) smoothly interpolates to the state of full equi- 
libration, Eq. ([?]) ■ when the length of the wire exceeds 
the equilibration length £ eq . The fully equilibrated dis- 
tribution Eq. ((7} is obtained from Eq. ^ by setting 
fi R = /j, L = fj, eq ; thus Afj, = 0, and also u = Vd- 

D. Brief summary 

The regime of partial equilibration described by the 
distribution function in Eq. ^ covers a wide range of 
lengths, £ a <C L < £ eq . It is more likely to be realized in 
experiments than the fully equilibrated regime ([7]) , as the 
length scale £ eq is exponentially large. Depending on the 
wire length L a particular state of the electron system 
is characterized by the extent to which the difference in 
chemical potentials Ap and temperatures AT of left and 
right movers has relaxed. The recent works Refs. I2IH22I 
addressed transport properties of wires with lengths in 
the range 

la < 4 < L ~ £ eq , (8) 

which covers the crossover from the partially equilibrated 
regime Eq. ([6]) to the fully equilibrated regime Eq. (J7J. 
The major emphasis of these works was on the effect of 
equilibration due to electron backscattering [Fig. HJd)]. 
The main focus of the present paper is on the transport 
properties of partially equilibrated wires with length 

e a 4ZL~ 4 « 4 9 . (9) 

In this regime the numbers N L and N R of the left- and 
right-moving electrons are individually conserved up to 



corrections as small as e~ M / T , so that electrons with en- 
ergies near the Fermi level pass through the wire with- 
out backscattering [Fig. [2c)]. This automatically implies 
that the conductance G unchaged intact by interactions 
and is still given by Eq. (|2|). However, electrons will expe- 
rience other multiple three-particle collisions [Fig. H{b)], 
which allow momentum and energy exchange within and 
between the two branches of the spectrum, thus altering 
thermal transport properties. The role of these processes 
was not explored in previous studies devoted to the trans- 
port in partially equilibrated quantum wires . 21 > 22 



III. BOLTZMANN EQUATION FORMALISM 

A. Three-particle collision integral 

Consider a quantum wire of length L, connected by 
ideal reflectionless contacts to noninteracting leads which 
are biased by a temperature difference AT (sec Fig. []}. 
In the following, we are interested only in the thermal 
transport properties of the wire and assume that there 
is no external voltage bias, V = 0. We describe weakly 
interacting one-dimensional electrons in the framework 
of the Boltzmann kinetic equation 

v p d x f(p,x)=l{f(p,x)}, (10) 

where v p — p/m is the electron velocity and the evolution 
of the distribution function is governed by the collision in- 
tegral X{f(j>,x)}. We consider the steady-state setup in 
which the distribution function does not depend explic- 
itly on time. The collision integral, in general, is a nonlin- 
ear functional of f(p, x), whose form is determined by the 
scattering processes affecting the distribution function. 
As discussed above, in our case the dominant processes 
are three-particle collisions. Assuming that the collision 
integral is local in space, we have 

P2.P3 
P 1 i,P 2 i,P 3 i 

x[/l/ 2 / 3 (l-/l0(l-/20(l-/3') 

-/i'/2'MW0(l -/*)(!- /a)], (11) 

where IF^f 3 is the scattering rate from the incoming 
states {pi,P2,P3} into the outgoing states {py ,p2< ,P3>}, 
and we used the shorthand notation = f(pi,x). The 
Boltzmann equation [Eq. (|10[) ] is supplemented by the 
boundary conditions stating that the distribution f(p, 0) 
of right-moving electrons (p > 0) at the left end of the 
wire and f(p, L) of left-moving electrons (p < 0) at the 
right end coincide with the distribution function in the 
leads, Eq. (g]). We note here that although Eq. (JTTJ) 
is written for the spinless case our subsequent analy- 
sis and solution of the Boltzmann equation presented in 
Sec. IIIIBllIIIDI is applicable to the spinful electrons as 
well. 
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An exact analytical solution of the Boltzmann equa- 
tion [Eq. ([TO])] is, in general, very difficult to find due 
to the nonlinearity of the collision integral Eq. (fTT|) . A 
simplification is, however, possible in the case of a linear 
response analysis in the externally applied perturbation 
(in our case, the temperature difference AT). Then the 
collision integral can be linearized near its unperturbed 
value. It is convenient to present f(p,x) as 



f(P,x) 



f 



(12) 



where f£ = (e i -- ep ~^/ T + 1) _1 is the equilibrium Fermi 
distribution function and ip(p,x) oc AT. When lineariz- 
ing Eq. (fTTj) with respect to tp(p, x) the factor /° (f — /°) 
in Eq. ([12")) makes it convenient to use the detailed bal- 
ance condition 

fpJlfpJ 1 - /£ )(i - /pJC 1 - O . ( 13 ) 



valid at e B 



Eq. (TT2"j) into the collision integral and using Eq. (TT3]) one 
arrives at the linearized version of Eq. (TIT]) 



T{^(pi,x)} = - 



P2 >P3 
P 1 /,P 2 /.P 3 / 



'l'2'3' 
H23 



X [i>(pi,x) +lj)(jP2,x) + 1p(p3,x) 

-ip(pv,x) - ip{p 2 ',x) ~ ip(p&,x)] ,(14) 



with the kernel K defined as 



^l'2'3' 
^123 



= W, 



1'2'3' 



123 



f0 rO fO 
J pi J P2-I P3 



*(i-0(i-/ P ° 2 ,)(i -/£,)• as) 



The explicit form of the scattering rate W{ 2 ^ is not 
important for the following discussion. It is discussed in 
detail in Appendix [B] 



B. Solution strategy 

Even after linearization the solution of the integral 
Boltzmann equation that satisfies given boundary con- 
ditions is still a complicated problem. However, our task 
is simplified greatly since we already know the structure 
of the distribution function f(p, x). Indeed, we have dis- 
cussed in Sec. [TT] that for wires with length £ a <C L ~ if, 
the electron system is in the regime of partial equili- 
bration with the distribution function given by Eq. ([3]). 
Thus, the class of functions we need to consider to solve 
our boundary problem is, in fact, rather narrow. The so- 
lution we are seeking is conveniently parametrized by six 
unknowns: [i R ( L '(x), u R ( L '(x), and T R ( L '(x). Instead 
of solving one intcgro-diffcrcntial equation for f(p,x) we 
will reduce our task to solving a system of six linear ordi- 
nary differential equations that govern the spatial evolu- 
tion of the parameters defining the distribution function 



in Eq. ([5]). This is possible since the momentum depen- 
dence of the distribution function is fully determined by 
our ansatz (|5|) , which allows completion of all p integra- 
tions in the Boltzmann equation analytically. Among the 
six equations we need, four represent conservation laws: 
conservation of numbers of right- and left-moving elec- 
trons N R / L , of total momentum P, and of energy E of 
the electron system. The other two are kinetic equations 
that account for the momentum and energy exchange be- 
tween subsystems of right and left movers, thus capturing 
the processes of thermalization. 

For the following analysis it is convenient to measure 
the chemical potentials and temperatures from their equi- 
librium values: n R{ - L \x) = ^+5^ R{L) {x) &ndT R< > L \x) = 
T + 8T R ^{x). Expanding now Eq. §5§ to the linear or- 
der in u R( - L \ Sn R(L) and 5T R{L) and using Eq. (|12[) we 
can identify ij;(p,x) that enters the collision integral in 
Eq. (JTID as 



ifj{p, x) = ip R (p, x) + V> L (p, x) 



£ p 2 > + £ p 3 >- Substituting wne re 



ip R M(j>,x) = d(±p) 

V£< £) (p, x) + ^ {L) (p, x) + < (L) (p, x) 

Here the three contributions are 

Sfi R ( L \x) 



T 

pu R ^ (x) 
T 



(16) 

(17) 

(18) 
(19) 



^,*)= (e '-y L)(g) . (20) 

These functions evolve in the real space as prescribed by 
the collision integral Eq. (|14|) while their boundary values 
can be extracted from the respective distributions in the 
leads [Eq. (gj]. Indeed, expanding Eq. (gj with V — 
one obtains 

f P = f° P + (£p ~r 2 )Ar / P °(i - OMp) 0{-p)\. (2i) 

By matching this result to Eq. (fT2|) with ip(p, x) taken 
from Eqs. ((TS]) and (JTTJ) we deduce the boundary condi- 
tions 



5n R (0) = 5[i L (L)=0, 

u R (0) = u L (L) = 0, 

5T R (0) = ST L (L) = AT/2. 



(22) 
(23) 
(24) 



Our task now is to derive the set of coupled ordinary 
differential equations that govern the spatial evolution 
of the unknown parameters u R ^(x), S (i R ( L ^ (x) , and 
ST R ( L ^(x). This will give us complete knowledge of the 
electron distribution function. Knowing all parameters 
in Eq. ([5]) we will be able to find the heat current and 
finally the thermal conductance of the system. Before we 
realize this plan, the conservation laws must be discussed. 
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C. Transport currents and conservation laws 



vation Eq. (|30j) . we find 



Conservation of the total number of particles implies 
that in a steady state the particle current I(x) is uni- 
form along the wire. Correspondingly, we infer from the 
conservation of total momentum P and total energy E 
that in the steady state a constant momentum current 
Ip and a constant energy current Ie flow through the 
system. In the following it will be convenient to express 
these currents as the sums of individual contributions of 
the left- and right- moving electrons, e.g., I = I R + I L , 
thus introducing: 

I R(L) (x) = | + °° f 0(±p)v p f(p, x) , (25) 
f +oc dv 

J f6(±p)pv p f(p,x), (26) 
f +oc dv 

J f6(±p)e p v p f(p,x). (27) 



I? L) (x) 



The positive sign in the step function corresponds to right 
movers and the negative one to left movers. Since we 
neglect small backscattering effects, the numbers of right- 
and left-moving electrons are conserved independently: 



N R / L = . 



(28) 



It follows then immediately from the continuity equations 
that particle currents are uniform along the wire, 



d x I R (x)=0, d x I L (x)=0. 



Similarly we present the conservation of total momentum 
and energy, 



d x [l?(x) + lji{x)] =0, 
d x [l§(x) + I^x)] =0. 



(30) 
(31) 



As the next step we express the currents (|2"5"|) - (j2~T)) in 
terms of the parameters defining the electron distribution 
function. Specifically, we use Eq. and Eqs. (|T7]) - (|2"0"]) 
together with the current definition in Eq. (|25[) and thus 
find from Eq. (|29|) 



d6n RiL) 



dx 



TPf- 



d u R ^ 
dx 



1 - 



7T 2 T 2 



7jT 4 T 4 



24^ 2 38V 



(32) 



When deriving this equation from Eq. ([25]) we had to 
carry out a Sommerfeld expansion up to the fourth order 
in T/fi <C 1- Note here that even though backscattering 
is neglected Sfj, R ^ must change in space to accommo- 
date the conservation laws for currents [Eqs. (l29l) - (f3"Tj) ]. 
and n R ( L \x) = \i + 5[i R ( L \x) coincide with ^ r n) = [i 
only at the ends of the wire [sec the boundary conditions 
Eq.' (EH)]. 

We then perform a similar calculation for the momen- 
tum and energy currents. From the momentum conser- 



T ( 5tt 2 T 2 
7n 2 T 2 \ fd5T R 



^ F 2/1 I 12 fl 2 J \ il.r dx 



du 1 



du l 



40 [i 2 ) \ dx 



dST L 
dx 



0, 



(33) 



while from the energy conservation, Eq. (|31l) . 



PF 



T 
2/t 



1 + 



Jllll\ ( duR dyL 
40 n 2 ) \ dx dx 

fdST R dST L 



\ dx 



dx 



+ 



= 0. 



(34) 



When deriving these two equations we also made use of 
Eq. (|32p to exclude the chemical potentials of the right 
and left movers. 



D. Scattering processes and kinetic equations 

Although the total momentum and energy currents Ip 
and Ie are conserved, such currents taken for the left 
and right movers separately, I^ L and Ig , are not. In- 
deed, the three-particle collisions shown in Fig. &b) in- 
duce momentum and energy exchange between counter- 
propagating electrons. Let us focus on a small segment 
of the wire between the positions x and x + Ax, where 



(29) < x < L. The difference Ip' (x + Ax 



Ip^ L (x) is 



equal to the rate of change of the momentum of right- 
moving electrons P R / L = p R / L Ax. Here p R l L is the rate 
per unit of length. As a result, the continuity equation 
for the momentum exchange reads 



d x [l R {x) - I L P (x)\ =2p 2 



(35) 



Here we used p L = —p R , which is ensured by the 
conservation of total momentum. In complete anal- 
ogy we can now relate the difference of energy currents 
I R {x+Ax) — I R {x) to the corresponding energy exchange 
rate E R = e R Ax, which gives us 



d x [l R (x)-I R (x)]=2eJ 



(36) 



where we also used e = — e guaranteed by the energy 
conservation. The right-hand side of Eqs. (j3"5)) and ([3^]) 
can be calculated from the collision integral of the Boltz- 
mann equation [Eq. (fT4")l ] . 

There are two basic processes which contribute to p R 
and e R . The first one includes two right movers that scat- 
ter off one left mover. This process is shown in Fig.[2Kd). 
The other process, when two left movers scatter off one 
right mover, is equally important. Keeping both terms, 
using Eqs. (fT4"]) and we find for the relaxation rates 
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(details of the derivation are given in Appendix [SJ 



and (|40p . (1411) can be reduced to the following form: 



-kF 
-kF 



Pf u 



R 



U 



2t 

fi5T 



v F 

R ST L 



T 

£ 



Vp(P3' - PS? 



Pl>0,p 2 >0,p 3 <0 
P,/ >0,p„/ >0,p„, <0 



-1'2'3 
H23 



(37) 

(38) 
(39) 



Having determined relaxation the rates we now return 
to Eqs. ([33)) and (|3f)|) . Computing the momentum and 
energy currents of right and left movers in the same way 
as we did in the previous section from Eqs. (|26p - (f2"?| and 
using the relaxation rates from Eqs. (l3"T|) - (|3"8"f we find 
two additional equations which describe the relaxation 
of momentum: 



ir 2 T 2 (du R du L 
F 12/i 2 \ dx dx 
7r 2 T / d8T R dST R 

6fi \ dx dx 
hk F u R - u L 



5tt 2 T 2 
7tt 2 T 2 

lo - "^ 2 " 



and energy 



TT 2 T 



Pf 



6^ 2 

TT 2 T f dST R 

3/i \ dx 



2 fdu R du L 
\ dx dx 
d8T L 



(40) 



dx 



J^T 2 

+ lo~^ 2 ~ 

2hk F ST R 



T 



(41) 



Equations fl32]h(|M|) and (0DJ)— (@TJ) , together with the 
boundary conditions Eqs. (l22l) - (f24)) . represent the 
closed system of six coupled differential equations 
whose solution fully determines the six parameters 
fi R / L , u R ' L , T R / L that define the electron distribution 
function ([5]). We now find these parameters explicitly. 
For that purpose let us introduce dimensionless variables 



u R ± u L 



Vp 



0± 



5T R ± 5T L 



T 



and the microscopic scattering length 



7r 3 r 4 

360/x' 



■ (vft) . 



(42) 



(43) 



Calculation of it, requires a detailed knowledge of the 
scattering rate, implicit in the kernel K} 2 g 3 , as a function 
of momenta transferred in a collision. In Appendix [B] we 
provide this information for the case of the three-particle 
collisions under consideration and in Appendix [C] find £b 
explicitly for the spinless and spinful cases. 

After some algebra the coupled equations (|3"3")h (|3~4")) 



86+ 


drj- 
dx 


dx 


dy+ 

dx 


dx 


<96L 


V- 


dx 


4 ' 


drj- 


0_ 


dx 





which contain two dimensionless parameters 
29V T 2 _Z^ 2 



1 



(44) 
(45) 
(46) 
(47) 



— — ■ ( 48 ) 



where corrections of higher order in T/fi <C 1 were ne- 
glected. The remaining two equations for the chemical 
potentials of right and left movers [Eq. (|32p] are not writ- 
ten here for brevity. The latter do not enter the heat 
current and thus are not explicitly needed. Equations 
P1)) - (|4T|) can now be easily solved, with the result 



2T 



a- + a+e L 



<{x) 



AT a+e x / e " + a-e^-^f 1 " 



2T 



a_ + a+e 



L/t b 



(50) 



_ AT fi_(e x / tb - 1) - p+( e ( L - x y ib - e L l lh ) 



v R (x) 



2T 



a- + a+e 



L/e„ 



(51) 



L/e b 



11 U!_ f U.+ K.~l 

(52) 

where a± = 1 ± a, /3± = 1 ± /? and 9 R / L = {9+ ± 0_)/2, 
^R/L _ ^ + _|_ rj_y2. This concludes our solution of the 
Boltzmann equation for three-particle collisions. 



IV. HEAT CURRENT AND THERMAL 
CONDUCTANCE 

Complete knowledge of the distribution function ([5]) 
allows us to compute physical observables. Specifically, 
we are interested in the thermal conductance K. For the 
latter we need to evaluate the heat current 



Iq{x) = Ie{x) - nl(x) 



(53) 



By using Eqs. ([5]) we carry out a Sommerfcld expan- 
sion for the particle and energy currents I and Ie from 
Eqs. ([25]) and ([27]). up to the fourth order in T/fj, < 1, 
and then find from the above definition [Eq. JSc 



n 2 T 2 
3h 



r)+{x) 



■9-{x) 



(54) 
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which is presented here in our notation defined in 
Eq. flU]). With the help of Eqs. QS])-^ it can be read- 
ily checked that Iq is uniform along the wire. This fact is 
a priori expected and follows from the conservation laws, 
which we already explored above. By knowing Iq we can 
finally find the thermal conductance K(L) — Iq/AT as 
a function of the wire length, 



K{L) tanh(L/24) + j3 



K a i 8tanh(X/24)+£ 



(55) 



This is the main result of our paper. Note here that Kq = 
ir 2 T/3h for the case of spinless electrons, whereas Kq is 
given by Eq. ([3]) for electrons with spin. The functional 
form of K(L) remains the same in both cases except for 
the expressions for 4, which we discuss below. Let us 
now analyze limiting cases of Eq. (f55"j) and discuss the 
microscopic form of the scattering length 4- 

Equation (|55[) interpolates smoothly between two dis- 
tinct limits. In short wires, L <C 4, from the expansion 
of Eq. (|55p one obtains for the interaction- induced correc- 
tion to thermal conductance, SK = K — Ko, the following 
result: 



5K(L) 
Ko 



30 /i 2 4 ' 



L < 



b ■ 



(56) 



In such short wires electrons propagate from one lead 
to the other, rarely experiencing three-particle collisions 
of the type shown in Fig. [2jb) . Thus their distribution 
function is approximately determined by that in the leads 
[Eq. 0]. Under this assumption one can adopt the strat- 
egy of Ref. [ljl applied previously for the calculation of 
conductance and thermopower in short wires, and treat 
the collision integral of the Boltzmann equation pertur- 
batively, thus neglecting effects of thermalization on the 
distribution function. Technically speaking, this corre- 
sponds to a lowest order iteration for the Boltzmann 
equation, which amounts to substituting distribution Q 
into the collision integral (ITi|) to calculate the correction 
to Iq. This perturbative procedure immediately repro- 
duces Eq. ([56"]). 

It is physically expected that in longer wires particle 
collisions Fig. 0(b) should have a much more dramatic ef- 
fect on the distribution function and thus thermal trans- 
port. Indeed, once full thermalization has been achieved 
for L ^> it, we find from Eq. ([55)1 that the correction to 
thermal conductance saturates: 



5K{L) 

Kn 



30 /i 2 



< L < 4 



(57) 



One interesting aspect of Eq. (|57|) is that SK is indepen- 
dent of the interaction strength. It means that no matter 
how weak the interactions are, for sufficiently long wires 
thermalization between right- and left-moving electrons 
is eventually established, which leads to saturation of SK. 
The behavior of SK(L) as a function of the wire length 
is summarized Fig. [3J 



0.006 - 



0.002 - 




- L 

5 2l b 



FIG. 3: (Color online) Interaction-induced correction to the 
thermal conductance of a clean quantum wire as a function of 
its length plotted for different values of temperature (from the 
bottom to the top curve): T/fi = 0.05, 0.1, 0.15, 0.2. For L < 
lb the correction scales with L and saturates to a constant 
value oc (T/fj,) 2 once 4 C i in accordance with Eqs. (|56|) 
and (T5T)) . 



The interaction strength, however, sets the length scale 
4 at which thermalization occurs. Its actual dependence 
on temperature is determined by the phase space avail- 
able for a three-particle collision to occur and by the de- 
pendence of the corresponding scattering amplitude on 
momenta transferred in a collision. For spinless electrons 
and Coulomb interaction we find (see Appendix [C] for the 
derivation and additional discussions) 



€j 1 ~ k F \i(k F w)(e 2 /nv F K) 4 (T/fif 



(58) 



where w is thw wire width and \i(z) = z 4 In 2 (1 / z) . In 
the case of spinful electrons, the scattering length changes 
to 



(- 1 ~ k F X 2 {k F w)(e 2 /hv F nf(T/y) ln 2 ( M /T) , (59) 

where X 2 (z) = ln 2 (l/z). Contrasting Eqs. flSSJ) and (|59"|) . 
one sees that the spin of the electron plays an important 
role since the inverse scattering length of spinful elec- 
trons is significantly larger, by a factor of (fi/T) 2 » 1. 
This is a manifestation of the Pauli exclusion principle. 
Indeed, for three-particle scattering to occur electrons 
must approach each other on a distance of the order of 
~ kp 1 . When electrons are spinless the Pauli exclusion 
suppresses the probability of such scattering. In contrast, 
the suppression is not as strong when the total spin of 
the three colliding particles is 1/2 since at least two elec- 
trons may have opposite spins while exclusion applies to 
the third particle. This technical point and the impor- 
tance of the exchange effect in the scattering amplitudes 
are discussed in more detail in Appendix [B] 
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V. DISCUSSION 

In this paper we studied the thermal transport prop- 
erties of one-dimensional electrons in quantum wires. In 
this system equilibration is strongly restricted by the 
phase space available for electron scattering and con- 
servation laws such that leading effects stem from the 
three-particle collisions. This is in sharp contrast with 
higher-dimensional systems where already pair collisions 
provide electronic relaxation. Although our theory is ap- 
plicable only in the weakly interacting limit, the results 
presented are still beyond the picture of the Luttinger liq- 
uid since three-particle collisions are not captured by the 
latter. We have elucidated the microscopic processes in- 
volved in electron thermalization and developed a scheme 
for solving the Boltzmann equation analytically within a 
linear response analysis. Our approach allows us to find 
the thermal conductance for arbitrary relations between 
the wire length and microscopic relaxation length [see 
Eq. (55])]. 

In order to establish a connection to previous wor k 21 i 22 
we emphasize that our solution of the kinetic equations 
and the result for thermal conductance presented in 
Eq. (|55[) rely on the simplifying assumption that electron 
backscattering can be neglected. This is a good approx- 
imation except for the case of very long wires, L > £ eq , 
where the small probability of backscattering ~ e _/i / T is 
compensated by the large phase space available for scat- 
tering to happen. Accounting for the backscattering pro- 
cesses, it was found in Ref. [2l| that for wires with length 



<-eq 



the thermal conductance is 



K{L) 



Ko 



L + L 



(60) 



This result gives only an exponentially small correction 
to the thermal conductance, 5K/K = —L/£ eq cx e~^/ T , 
in the limit L <C £ e q, since in the analysis of Ref. [2l| 
thermalization effects on the distribution function were 
neglected. It is our result Eq. ((57)) that gives the leading- 
order correction to 5K in this case. On the other hand, 
our expression (|55[) is not applicable for the long wires, 
L ~ £ eq , whereas Eq. (f6"U|) works in this regime. It dis- 
plays an essentially different feature, which is solely due 
to backscattering processes, namely, the vanishing ther- 
mal conductance 5K oc 1/L as L — ¥ oo. 

Our work may be relevant for a number of recent 
experiments. In particular, the authors of Ref. [l^ re- 
ported thermal conductance measurements and a lower 
value of K than that predicted by the Wiedemann- 
Franz law, at the plateau of electrical conductance. As 
we explained, corrections to G are exponentially small, 
G = 2e 2 /h - 0(e-^ T ), for wires with L < £ eq . Thus 
the conductance remains essentially unaffected by inter- 
actions, and its quantization is robust. In contrast, the 
effect of three-particle collisions on the thermal conduc- 
tance is much more pronounced. Our equation (JSJJ shows 
that the thermal conductance is reduced by interactions, 



which is qualitatively consistent with the experimental 
observation^ Apparent violation of the Wiedemann- 
Franz law is due to the fact that interaction-induced 
corrections SK and SG originate from physically dis- 
tinct scattering processes [see Figs. [2jb) and Hfc), re- 
spectively] . 

Another experiment 1 ^ reported measurements of the 
electron distribution function in one-dimensional wires. 
This experiment demonstrated that electrons thermalize 
despite the severe constraints imposed by the conserva- 
tion laws and dimensionality on the particle collisions. 
We take the point of view that three-particle collisions 
are responsible for relaxation and provide an explicit so- 
lution of the Boltzmann equation, thus uncovering the 
structure of the distribution function [see Eqs. ([5]) and 
(|49l) - (f52|) ]. which in principle can be compared to exper- 
imental results^ 

A related studyi£ provided us information about the 
time scales of thermalization of one-dimensional elec- 
trons. Although we do not study the latter our results for 
the relaxation lengths Eqs. (|58p and (|59p can be directly 
linked to the experiment. Note also that the dramatic 
difference between the relaxation lengths, and thus the 
times, of spinful and spinless electrons provides a distinct 
signature of three-particle collisions that could be tested 
experimentally. 

There is a very important limitation on the applica- 
bility of Eqs. ([5"5f and ([ST?]) that we need to discuss in 
the case of spinful electrons. 24 From the point of view 
of Luttinger liquid theory, electrons are not well-defined 
excitations in one dimension and instead one should use 
a bosonic description in terms of charge and spin modes. 
The weakly interacting limit considered here and usage of 
the Boltzmann equation assumes that electrons maintain 
their integrity during collisions and thus neglects effects 
of spin-charge separation. In order to quantify to what 
extent such a description is valid, consider an electron 
with excitation energy £ above the Fermi energy fi. For 
quadratic dispersion, e p = p 2 /2m, the velocity of such 
electrons differs from that of the electrons in the Fermi 
sea by an amount Av — t^/mvp. Spin and charge do not 
separate appreciably if Av 3> v c — v s , where w c r s ) are 
the velocities of charge (spin) excitations. At finite tem- 
peratures the characteristic excitation energy is £ ~ T, 
so that the above condition can be equivalently reformu- 
lated as T / ji 3> (v c — v s )/vf- For weakly interacting 
electrons the difference between the velocities of charge 
and spin modes is related to the zero-momentum Fourier 
component of the electron-electron interaction potential, 
namely, v c — v s ~ Vq/ttK 1. This implies that at low 
temperatures when T//i <C Vo/fovp a description in terms 
of electrons breaks down and Eqs. (|55|) and (|59)l are no 
longer applicable. 

Finally, our work also points to open issues and direc- 
tions for future research. It is of great interest to un- 
derstand the fate of energy relaxation and the nature of 
thermal transport in the case of strong interactions which 
simultaneously have to be combined with nonequilibrium 
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conditions. At very low temperatures a description of a 
one-dimensional system in terms of electronic excitations 
becomes inadequate even if the interactions are weak. 
The effect of spin-charge separation has to be included, 
and thermal transport from plasmons and their relax- 
ation are central issues to consider. 
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Appendix A: Derivation of P R and E R 



In this appendix we derive Eqs. (|3"7| - (l3"9l presented in 
the main text of the paper. As explained in Sec. IIII D[ 
when computing P R and E R we have to account for two 
types of scattering process. One is shown in Fig.^b) and 
the other is similar and consists of a scattering of one left 
mover and two right movers. We start by considering the 
quantity V n — Y^ Pl >aPifl- Its rate of change is 



a, = 5>?a = - £ 

Pi>0 



Pi ^123 



PI >0,P2 1P3 
P 1 /,P 2 /.P 3 / 



X (ipi +Ip2 + lp3- IpV ~ ifo' ~ 1p3>) , 



(Al) 



where we used the Boltzmann equation [Eq. (fTU]) ] and 
the short-hand notation ipi — ip(pi,x). It is convenient 
to split each sum from the last equation into parts that 
contain positive and negative values of the momenta, so 
that one gets 

v n = £(...) + 3 £(...) + e £(...) + £(•••) 

+-- +-- ++- 

++- +-- +++ 

+ £(... ) + 3£(...) + 3£(...) + 2£(...) 

+++ +++ +++ ++- 

+-- ++- +++ 

+2£(...)+3£(...) + 6£(.. .)+£(■■•)■ 



++- 
++- 



+++ 
+++ 



The notations here are as follows 

£(-) = £ (•••) 



(A2) 



(A3) 



P1>0,P2<0.P3<° 
p 1 / <0,p 2 / <0,p 3 / <0 



and analogously for the other terms. When deriving 
Eq. (|A2[) we have used the following symmetry proper- 
ties of the kernel: (a) exchange of incoming and out- 



going momenta K{ 23 



1'2'3' _ w-123 



1'2'3'j O 3 ) pairwise exchange 



^123 3 — ^213 3 i an d ( c ) inversion of momenta p, — > —pi, 
K} 2 | 3 = Kljj~ 2 2 _3 3 . In the final expression for V n we 
keep only the terms of Eq. (| A2|) that contain equal num- 
bers of positive incoming and outgoing momenta, i.e., 
the last three terms. The other terms contain at least 
one state near the bottom of the band and therefore give 
a contribution that is exponentially suppressed due to 
the small probability of finding an unoccupied state. Af- 
ter employing Eqs. (|16p - (|20[) combined with momentum 
and energy conservations, we end up with 



V n = 3 £ K# 3 ' 3 ' 



++- 
++- 



ST R - ST 1 



2^2 



[2p" - (-P 3 ) n ] 



x(e p 



T 



[2p? + (-p 3 ) n ](P3-P3') 



For n = 1 from the last expression we easily get 

L E^f'fe'-K) 2 , (A4) 



P r = Vi = -3- 



T 



++- 

++- 



which reduces to Eq. (|37p in the main text. For n — 2 



2m T 2 



£ K 123 3 ( e P 3 - ~ £ p 3 ) 2 



++- 
++- 



(A5) 



One additional step is required to obtain Eq. (|3"5)) . Since 
all three particles participating in a collision are located 
near the Fermi points it means that the characteristic 
momentum of right movers is ~ pf while for the left 
mover it is ~ —pf- In contrast, the momenta transferred 
in a collision qi = py — pi are much smaller ~ T/vf <C 
Pf, which stems from the temperature smearing of the 
occupation functions implicit in the kernel K} 2 3 3 ■ Since 
\Qi\ Pf we approximate p^ w — pp and linearize the 
spectrum near the Fermi points, in particular, 



1 

2m 



[(-PF + q 3 ) 2 



Pf] 



-v F q 3 , (A6) 



which then brings the last expression for E to the form 
of Eq. (|38p in the main text. 



Appendix B: Three-particle scattering amplitude 

For most of our analysis the detailed form of the scat- 
tering rate entering kinetic equation (|10|) was not im- 
portant. However, for the calculation of the microscopic 
quantities, such as the scattering lengths £ a and lb, we 
need to know the precise form of the scattering rate 
W123 3 introduced in Eq. (Tn|) . Below we give the de- 
tails of the structure of the scattering rate. We start 
with the golden rule expression 



W-, 



1'2'3' 



2tt 



123 



\A\ 2 



2 5{E - E')6 p ,p> 



(Bl) 



11 



a) 



A(ll',22',33') 

W Pr 



b) 



A(12',21',33') 




FIG. 4: (Color online) Direct (a) and five exchange (b)-(f) 
terms in the three-particle amplitude A^23 3 [Eq. (|B2|) ] that 
contribute to the finite momentum P R and energy E R ex- 
change rates between right and left movers. 



where A\ 2 \ 3 is the corresponding scattering amplitude, 
while the 5-functions impose conservations of the to- 
tal energy E(E') = 2~2i £ Pi(p-i) anc * total momentum 
P(P') = J2iPi(i') °f the colliding electrons. One should 
note that in Eq. (|B1[) we include 6p t p> in the definition 



of the scattering rate W^i 3 ra ther than the amplitude 
A123 3 , which is in contrast to the usual convention. This 
step simplifies our notations. 

Since the electrons interact with two-particle interac- 
tion potential V{x), the three-particle scattering ampli- 
tude is found to the second order in V(xY The details of 
this calculation were presented in Ref. |l9|. In the case of 
spinless electrons the final result reads 



A 



1'2'3' 
123 



sgn(l'2'3')A(ll',22',33'). (B2) 



tt(1'2'3') 



One should notice that Eq. (|B2I) contains the term 
A(lV, 22', 33'), which is the amplitude of the direct scat- 
tering process [Fig.UJa)], and the terms obtained by the 
remaining five permutations of the outgoing momenta, 
which are the exchange terms [Figs. IUJb)-[2ff)]. They 
can be written compactly for the segment Ax of the wire 



as follows: 



A(la,2b,3c) = af 2 + < 3 C + a b 2 c 3 , 



(B3) 



af 7 = at"™ 



1 



7 — v , Vrt, 



~ (Air) 2 Pa ~ Pl pi>-p 2 
1 



E £ pi £ pb Ep—p 1 -p l 
1 



E—e Pa —e P2 —sp- 



Pa~P2 



, (B4) 



where (a, b, c) is a particular permutation of (1', 2', 3'). In 
Eq. (|B2[) the notations 7t(. . .) and sgn(. . .) denote permu- 
tations of the final momenta and parities of a particular 
permutation. Finally, V p is the Fourier-transformed com- 
ponent of the bare two-body interaction potential. For 
the calculations we take the Coulomb interaction between 
electrons, 



V(x) = - 



\Jx 2 + 4w 2 \Jx 2 + Ad 2 



(B5) 



screened by a nearby gate, which we model by a conduct- 
ing plane at a distance d from the wire. We also intro- 
duced a small width w of the quantum wire, w <C d, to 
regularize the diverging short-range behavior of this po- 
tential. This enables us to evaluate the small-momentum 
Fourier components V p of the interaction potential V(x). 
To this end, we find in the limit h/d -C p <C h/w 



Vn 



2e 2 



In 



Pw 

\p\ 



1 



P 



while in the limit of very small momenta p -C h/d 



K \W 



p 2 ln(p d /\p\) 



p 2 ]n(d/w) 



(B6) 



(B7) 



In the last two equations we introduced the notations 
Pw = h/w and pd = h/d. We also employed logarith- 
mic accuracy approximation for V p , meaning that nu- 
merical coefficients in the arguments of the logarithms in 
Eqs. (|B6[) and (|B7[) are neglected. In the following dis- 
cussions we refer to Eq. (|B6|) as the unscreened Coulomb 
potential and to Eq. (|B7|) as the screened one. The com- 
plete expression for the amplitude (IB2[) with the inter- 
action potential taken in the form (|B6I) or (|B7|) is fairly 
complicated. However, major simplification is possible by 
studying the kinematics of the three-particle collisions, 
which in a way allows us to obtain the approximated 
form of the amplitude for specific scattering processes, 
such as the one in Fig. EJb), which determines the scale 
Ih- 

It is convenient to label the outgoing momenta as 
Pi' = Pi + Qi for i = 1, 2, 3 in order to separate explicitly 
the momenta qi transferred in a collision. Momentum 
conservation then reads 



qi + 92 + 93 = . 



(B8) 
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while energy conservation E 
rewritten as 

Ipiqi + 2p 2 q 2 + 2p 3 q 3 



E' can be equivalently 



4 



o. 



(B9) 



At low temperatures, T -C /i, the Fermi occupation func- 
tions constrain particles participating in the collision to 
lie in a momentum strip of the order of T/vp <C pp near 
the Fermi level. This means in practice that the typi- 
cal momentum transferred in a collision will not exceed 
max{|g|} < T/vp. To leading order in T/p <C 1 the en- 
ergy and momentum conservation requirements for the 
scattering process in Fig. 0] can be resolved by 



qi « ~qi + 0{[(p! - p 2 ), q2\/PF} 



(BlOa) 



and 



93 



gi(gi +Pi -P2) 
2p F 



0{[( Pl ~p 2 ), qi } 2 /p 2 p}, (BlOb) 



where we used p\ — p 2 ~ T jvp and set pz ~ ~Pf- From 
this analysis one concludes that energy transfer between 
the right and left movers occurs via small portions of 
momentum q 3 exchange such that 

{M, \<h\} ~ T/vp , |g 3 | ~ T 2 /v Ft i « {\ qi \, \q 2 \} . 

(Bll) 

Having two small parameters at hand, *C 1 and 

1 93 1 /1 1i I ^ 1) an d accounting for all the exchange con- 
tributions, one can expand the amplitude (|B2j) to the 
leading nonvanishing order^ In the course of this ex- 
pansion we observed that exchange contributions result 
in severe cancellations between different scattering pro- 
cesses. The result of the calculations for the model of 
unscreened interaction (IB6I) is 



2e^ 



9Ai(/cfw) 
64^ 2 (Aa;) 4 



(B12) 



There are several general remarks we need to make re- 
garding the scattering amplitude in Eq. (|B2I) . It is known 
from the context of integrable quantum many-body prob- 
lems^ that for some two-body potentials, TV-body scat- 
tering processes factorize into a sequence of two-body col- 
lisions. In the context of this work, this means that three- 
particle scattering for the integrable potentials may result 
only in permutations within the group of three momenta 
of the colliding particles; all other three-particle scatter- 
ing amplitudes must be exactly zero for such potentials. 
We have checked explicitly that the three-particle scatter- 
ing amplitude in Eq. (|B2[) is nullified for several special 
potentials: for the contact interaction, V p — const, for 
the Calogero-Suthreland model, V p oc \p\, and also for the 
potential V p oc \—p 2 jp\ which is dual to the bosonic Lieb- 
Liniger model. Surprisingly, we have also noticed that 
the logarithmic interaction potential VJ, oc In \p\ gives ex- 
actly zero for the three-particle amplitude in Eq. (|B2[) 
although we are unaware of any exactly solvable model 
for that case. This is the reason to keep the next leading- 
order term ~ (p/Pw) 2 <C 1 in Eq. (1B6[) . which prevents 
the amplitude in Eq. (|B2I) from vanishing exactly. 

The second set of remarks concern electrons with spin. 
For the latter the three-particle amplitude has the same 
form as Eq. (|B2|) . but it acquires an additional depen- 
dence on the spin indices: 



= E sign(l'2'3')S^ CT3 ,A(ll',22',33'), 



7r(l'2'3') 



(B14) 
One can repeat 



where -"i^^ — A A A 
the expansion of the amplitude for \qi\/pF <C 1 and 
l^sl/lffil ^ 1 an d observe that due to the spin structure 
the exchange terms do not cancel each other. In partic- 
ular, with the help of Eq. (|B6I) we find the amplitude for 
the case of the unscreened Coulomb potential in the form 



where the amplitude is written for a segment of the wire 
of length Ax and the function \i(kpw) was introduced 
earlier [see the definition after Eq. (|58|) ]. For the screened 
case we find 
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2e 2 \ 25A 3 (fc F ef) 



K J V(Ax) 4 



4 In 

Pf 



PF 



Ill 
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, (B13) 



where A 3 (z) = z 8 ln 2 (l/z). Both amplitudes (|B12|) and 
(IB13|) are written in logarithmic accuracy approximation. 
As argued above, the typical scattering processes studied 
here involve only small-momentum transfer, of the order 
of q ~ T/vp. As a result, from the conditions of appli- 
cability of the interaction potential Eq. (|B6I) it follows 
that the corresponding amplitude Eq. (|B12I) applies for 
T 3> hvp I d. Similarly, the screened interaction potential 
Eq. (|B7I) and corresponding amplitude Eq. (IB13|) apply 
at lower temperatures T hvpjd. 
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2.1-4123 I -{—) 3 2 M 2( Aa; )4 
W 

(B15) 

which is by a factor of (pf/|qi|) 2 1 larger than 
Eq. (|B12[) : \ 2 (k F w) was defined under Eq. (1531) . 



Appendix C: Intra-branch and inter- branch 
relaxation lengths 



In this appendix we estimate the scattering lengths £ a 
and £b- Our starting point for evaluation of the inter- 
branch length £b is the expression 



1080m 4 



1 



7T 3 T 4 VFkpAx 



E 

++- 

++- 



(vf qzf 



ir, 



1'2'3' 
123 



n/ }, 



(Cl) 



13 



which follows from Eqs. (|39[) and (|43[) . where in addition 
we introduced the notation 

^{/°} = — /pi+gJ/paO' — fp2+q2)fp3^ ~ fp3+Q3^- 

(C2) 

In view of the kinematic constraints (jBlOal) and (jBlObl) , 
conservation of momentum and energy in the expression 
(IB1[) for the scattering rate 3 can be presented as 



6{E-E')5p, P , « - — S \ q 3 
2vp 



gl(gl +Pl -P2; 

2p F 



(C3) 

which eliminates two out of six momentum integrations 
in Eq. (|C1|) . The other four integrals can be completed 
analytically with logarithmic accuracy. This amounts 
to replacing the weak logarithmic parts of the ampli- 
tude in Eqs. (lB12|) . (jBTl . and (|BTo|) by their typi- 
cal values taken at characteristic momenta q\ ~ T/vp 
and ~ T 2 jvp\i. We thus treat ln(q 2 /2pF|g 3 |) in 
Eq. (|B12|) as a constant of order unity and approximate 
ln(pf/|«i|) ~ ln(|gi|/|g 3 |) ^ ln(/*/T) in Eqs. flEE^ and 
(|B15j) . After this step we can integrate in Eq. (|C1|) explic- 
itly by linearizing the electron dispersion relation inside 
the Fermi functions and get 



E" 4 
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(C5) 



P1P2P3 \ 21 / 

For the spinless case and high-temperature regime T 3> 
Hvp/d, where the Coulomb interaction is unscreened, we 
obtain 



Xi(kpw) ( e 2 \ 
p 2 F T A Ax 



hn J 



(4 



sinh 



2 I vpqx ' 
2T , 



(C7) 



Note here that we do not keep track of the numerical co- 
efficient in the expression for £{, since within the adopted 
calculation with logarithmic accuracy this coefficient is 
not determined. After the remaining qi integration one 
recovers Eq. (j58")l . presented in the main text of the paper. 

At lower temperatures T <C hvp/d, screening effects 
become important and one should use Eq. (|B13p in the 
expression for the scattering length Eq. (|Cip . Estimation 
of £b in this case gives 



r 1 



k F X 3 (kpd)(e 2 /hv F K)\T/ a*) 7 In 2 (p/T) . (C8) 



In the spinful case this calculation is completely anal- 
ogous to the one above; we just need to use a different 
expression for the scattering amplitude. With the help 
of Eq. (|B15[) . which is applicable for the model of an 
unscreened Coulomb potential interaction, we get at the 
intermediate step with logarithmic accuracy, 



! _ \ 2 {k F w) t <?\ 2 / £ \ ^ -a y 2 ' vi 
b ~ T^Ax \hKj m \T)^ sinh 2 (^) 



(C9) 



After the final integration this translates into Eq. (|59|) . 

We turn now to discussion of the intra-branch relax- 
ation length £ a introduced in Sec. [TTJ Unlike the case of 
interbranch relaxation, Fig. [TJb), here all three colliding 
particles are near the same Fermi point; see Fig. [Ha). 
In this case, the typical momentum change for the three 
electrons is the same, 



ml 



\q 2 \ ~ \q 3 \ ~ T/vp. 



(CIO) 



At this point we should emphasize that for the processes 
that determine the length scale £b, a new energy scale 
T 2 //i appeared in the problem purely from the kinematic 
constraints based on the conservation laws. This scale 
determined the typical momentum transfer of the parti- 
cle that was alone at one side of the Fermi surface; see 
Eqs. (|B10bp and [[BlT]) . 

Another important quantity is the scattering ampli- 
tude. For Coulomb interaction and for the process where 
all three particles are near the same Fermi point, it 
is a relatively complicated expression, but similarly to 
Eqs. (fB~T2|) and (|B13|) it depends on momenta only weakly 
(logarithmically) for the intra-branch processes. 

These two observations help us to estimate £ a using 
the known result Eq. (|58p for l^. Namely, by replacing 
the energy scale T 2 / \i in £\, by T, we obtain the estimate 



kp\i(kpw)(e 2 /HvpK,) 4 (T/iJ,y 



(Cll) 



for the unscreened Coulomb case, T 3> Hvp/d. At lower 
temperatures, T <§C hvp/d it changes to £~ x cx T 6 . It is 
important to emphasize that regardless of the interaction 
model we use there exists a distinct separation between 
the scales of the relaxation lengths, namely, 



la/tb ~ T/fX « 1. 



(C12) 



This fact justifies our ansatz for the distribution function 
[see the discussion after Eq. ([SJ]. The detailed calculation 
of £ a will be presented elsewhere^ 



1 A. A. Abrikosov, Fundamentals of the Theory of Metals 2 G. V. Chester and A. Thellung, Proc. Phys. Soc. London 
(North-Holland, Amsterdam, 1988). 77, 1005 (1961). 



14 



3 Electron- Electron Interaction In Disordered Systems, 
edited by A. J. Efros and M. Pollak (Elsevier, Amsterdam, 
1985). 

4 G. Catelani and I. L. Aleiner, JETP 100, 331 (2005). 

5 T. Giamarchi, Quantum Physics in One Dimension, 
(Claredon Press, Oxford, 2003). 

6 D. L. Maslov and M. Stone, Phys. Rev. B 52, R5539 
(1995); V. V. Ponomarenko, Phys. Rev. B 52, R8666 
(1995); I. San and H. J. Schulz, Phys. Rev. B 52, R17040 
(1995). 

7 R. Fazio, F. W. J. Hekking, and D. E. Khmelnitskii, Phys. 
Rev. Lett. 80, 5611 (1998). 

8 C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68 1220 
(1992); 76, 3192 (1996). 

9 K. Schwab, E. A. Henriksen, J. M. Worlock, and 
M. L. Roukes, Nature Phys. 404, 974 (2000). 

10 O. Chiatti, J. T. Nicholls, Y. Y. Proskuryakov, N. Lump- 
kin, I. Farrer, and D. A. Ritchie, Phys. Rev. Lett. 97, 
056601 (2006). 

11 J. T. Nicholls and O. Chiatti, J. Phys.: Condens. Matter 
20, 164210 (2008). 

12 F. Sfigakis, A. C. Graham, K. J. Thomas, M. Pepper, 

C. J. B. Ford, and D. A Ritchie, J. Phys.: Condens. Matter 
20, 164213 (2008). 

13 G. Granger, J. P. Eisenstein, and J. L. Reno, Phys. Rev. 
Lett. 102, 086803 (2009). 

14 C. Altimiras, H. le Sueur, U. Gennser, A. Cavanna, 

D. Mailly, and F. Pierre, Nature Phys. 6, 34 (2009). 

15 Y. F. Chen, T. Dirks, G. Al-Zoubi, N. O. Birge, and N. Ma- 
son, Phys. Rev. Lett. 102, 036804 (2009). 

16 G. Barak, H. Steinberg, L. N. Pfeiffer, K. W. West, 
L. I. Glazman, F. von Oppen, and A. Yakoby, Nature Phys. 
6, 489 (2010). 

17 D. B. Gutman, Y. Gefen, and A. D. Mirlin, Phys. Rev. 
Lett. 101, 126802 (2008); Phys. Rev. B 81, 085436 (2010). 



18 S. Takei, M. Milletari, and B. Rosenow, Phys. Rev. B 82, 
041306 (2010). 

19 A. M. Lunde, K. Flensberg, and L. I. Glazman, Phys. Rev. 
B 75, 245418 (2007). 

20 J. Rech, T. Micklitz, and K. A. Matveev, Phys. Rev. Lett. 
102, 116402 (2009). 

21 T. Micklitz, J. Rech, and K. A. Matveev, Phys. Rev. B 81, 
115313 (2010). 

22 A. Levchenko, T. Micklitz, J. Rech, and K. A. Matveev, 
Phys. Rev. B 82, 115413 (2010). 

23 K. A. Matveev, A. V. Andreev, and M. Pustilnik, Phys. 
Rev. Lett. 105, 046401 (2010). 

24 T. Karzig, L. I. Glazman, and F. von Oppen, Phys. Rev. 
Lett. 105, 226407 (2010). 

25 A. Levchenko, Z. Ristivojevic, and T. Micklitz, Phys. Rev. 
B 83, 041303(R) (2011).' 

26 Y. M. Sirenko, V. Mitin, and P. Vasilopoulos, Phys. Rev. 
B 50, 4631 (1994). 

27 To be accurate we point out that for the noninteracting 
electrons Go = (2e 2 //i)(l + e _My/T ) — ^ where the exponen- 
tially small term comes from the states at the bottom of the 
band. At low temperatures e _M//T <C 1 and thus all such 
small contributions are neglected throughout the paper. 

28 For the estimates we use typical parameters from the ex- 
periment of Ref. HH . In particular, tip ~ 8 x 10 5 m/s, 
kp ~ 10 8 m _1 , k ~ 10, and w ~ 5 nm. This translates 
into e 2 /HvfK ~ 1 and fj, = /wf&f/2 ~ 25 meV so that for 
the temperature T ~ 1 K the thermal relaxation length 
[Eq. ((SSJl] is of the order l b ~ 1 /im. 

29 Z. Ristivojevic, A. Levchenko and K. Matveev, (unpub- 
lished) . 

30 B. Sutherland, Beautiful Models (World Scientific, Singa- 
pore, 2004), Sees. 1.3-1.5. 



